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The relations between matrix elements of different P- and T-odd weak interactions are derived. 
We demonstrate that similar relations hold for parity nonconserving (PNC) transition amplitudes 
and electron electric dipole moments (EDM) of atoms and molecules. This allows to express P- and 
T-odd effects in many-electron systems caused by different symmetry-breaking mechanisms via each 
other using simple analytical formulas. We use these relations for the interpretation of the anapole 
moment measurements in cesium and thallium and for the analysis of the relative contributions 
of the scalar-pseudoscalar CP-odd weak interaction and electron EDM to the EDM of Cs, Tl, Fr 
and other atoms and many polar molecules (YbF, PbO, ThO, etc.). Model-independent limits on 
electron EDM and the parameter of the scalar-pseudoscalar CP-odd interaction are found from the 
analysis of the EDM measurements for Tl and YbF. 
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I. INTRODUCTION 

The study of the parity and time invariance viola- 
tion in atoms, molecules and nuclei is a low-energy, rel- 
atively inexpensive alternative to the high-energy search 
for new physics beyond the standard model (see, e.g. a 
review [i| ). Accurate measurements of the parity non- 
conservation (PNC) in atoms is one of the most promising 
ways of exploring this path. It culminated in very precise 
measurements of the PNC in cesium Q ■ Interpretation of 
the measurements based on accurate atomic calculations 
led to perfect agreement with the standard model and 
put strong constrains on any new physics beyond it 0, H| 
(see also review [1| for more detailed discussion). First 
unambiguous measurement of the nuclear P-odd anapole 
moment was also reported in the Cs PNC experiment Q ■ 

At present, the study of symmetry violations in atoms 
and molecules goes mostly in three major directions (see, 
e.g. [fl): (i) the PNC measurements for a chain of iso- 
topes; (ii) the measurements of nuclear anapole moments; 
and (iii) the measurements of the P,T-odd permanent 
electric dipole moments of atoms and molecules. Inter- 
pretation of the anapole moment and EDM measure- 
ments requires sophisticated atomic or molecular cal- 
culations. The calculations are difficult and sometimes 
disagree with each other. For example, the calcula- 
tions of the thallium electric dipole moment (EDM) by 
Nataraj et al [(| and by Sahoo et al Q disagree with ear- 
lier calculations [1, Q , calculations of the nuclear spin- 
dependent PNC in cesium by Mani and Angom [l(| dis- 
agree with earlier calculations pd| - [T3 |. etc. The diffi- 
culties are mostly due to the strong inter-electron cor- 
relations which need to be treated to all orders of the 
many-body perturbation theory. Majority of publica- 
tions present the effect of one particular P- or CP-odd 
interaction for a specific atom or a molecule. If another 
weak interaction Hamiltonian is considered, the calcula- 
tions are to be done again. 



The extraction of the nuclear anapole moment from 
the PNC measurements for thallium [l5( also presents a 
problem. The analysis of the experimental data based 
on atomic calculations favors a negative value for the 
thallium anapole moment [TBI . ITfjl w hile all nuclear cal- 
culations produce positive values [171 ]. 

Since the analysis based on sophisticated atomic cal- 
culations is very complicated it is important to have ad- 
ditional tools which may help to check the calculations 
for consistency or even avoid the calculations. In present 
paper we build such tools by studying the relations be- 
tween matrix elements for different symmetry-breaking 
operators. We show that using these relations the result 
for a specific operator can be easily obtained analytically 
if the calculations for another operator are available. The 
relations can also be used to check different calculations 
for consistency. Additional benefit comes from the need 
to separate the effects of different sources to the PNC or 
EDM of atoms or molecules in the analysis of the exper- 
imental data. Recent measurement of T and P violation 
in YbF molecule [l8[ combined with the Tl EDM mea- 
surement [l9j would allow to obtain independent limits 
on electron EDM and CP-violating interaction. However, 
the relative sign of these contributions is different in dif- 
ferent calculations [H, [2(| ■ This sign problem is solved in 
the present work by calculating the ratio of the matrix 
elements. Then independent limits on electron EDM and 
scalar-pseudoscalar CP-odd interaction is presented. 

There is an active search for the CP-violating per- 
manent EDM of polar molecules, such as YbF [rsf . 
PbO [H [H[, ThO [H, etc. Interpretation of the 



measurements requires complicated molecular calcula- 
tions. Most of calculations consider only one possible 
source of the molecular EDM: electron EDM or scalar- 
pseudoscalar CP-odd interaction. We show that the ana- 
lytical ratios of the matrix elements provide reliable link 
between molecular EDMs caused by different sources. 
Therefore, if calculations exist for a particular CP-odd 
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operator, no new calculations are needed to find an EDM 
caused by a different operator. 

The approach developed in this paper is easy to apply 
when matrix elements of the symmetry-breaking opera- 
tor between s l / 2 and p 1 / 2 single-electron states strongly 
dominate over other matrix elements. This is the case for 
the PNC amplitudes between atomic s-states and for the 
EDM of atoms and molecules caused by electron EDM 
or scalar-pseudoscalar CP-odd weak interaction. We also 
use this approach for the thallium anapole sign problem. 



II. SINGLE-ELECTRON MATRIX ELEMENTS 
OF WEAK INTERACTION AND MANY-BODY 
EFFECTS 

The PNC amplitude of an electric dipole transition 
between states of the same parity \a) and \b) is equal to: 
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where d = — e J2i r i is the electric dipole operator, -ffpNC 
is the operator of a P-odd CP-even weak interaction. 
The EDM of an atom in state a is given by 
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where Hep is the operator of a CP-odd weak interaction. 

States a, b and n in ((T|) and @ are the many-electron 
atomic states. However, we can start for simplicity from 
an atom with one external electron above closed shells 
keeping in mind cesium atom as an example. Then in 
the lowest order of the perturbation theory in residual 
Coulomb interaction the PNC amplitude and the EDM 
of the atom are given by (TTJ) and ^ in which states a, 
b and n are single electron (e.g. Hartree-Fock) states. It 
is easy to see that all single-electron matrix elements of 
weak interaction are proportional to each other. Indeed, 
only short distances contribute to the value of the ma- 
trix elements (see appendix for analytical estimations). 
These distances are r < for the PNC interactions 
and for the scalar-pseudoscalar CP-odd interaction (Rn 
is nuclear radius) and r < ao/Z for the electron EDM op- 
erator (ao is Bohr radius). On these distances the energy 
of the single-electron state can be neglected compared to 
the nuclear potential —Ze 2 jr and Dirac equations for all 
single-electron states become identical. The only differ- 
ence comes from the normalization of the states. Suppose 
we have two different operators of, say, CP-odd weak in- 
teraction i?cp and -ffcp- If we establish proportionality 
between matrix elements for a particular pair i, j of single 
electron states 



TABLE I: The ratio of the matrix elements of the spin- 
dependent to spin- independent P-odd weak interactions (i?i) 
and scalar-pseudoscalar CP-odd interaction to electron EDM 
(-R2) for the 6s,6pi/2 states of cesium with and without the 
inclusion of dominating many-body effects. The ratios are 
stable while the matrix elements change by up to two times. 
Units: x/(-Q w ) for Ri and 10~ 15 C PS /d e a.u. for R 2 . 
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"Brueckner orbitals 
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then due to the proportionality of the wave functions 
on short distances, the same proportionality, with the 
same value of R would hold for any pair of single-electron 
states and for the total EDM of atom ([2]). Furthermore, 
if the S1/2 — P1/2 matrix elements of the weak interaction 
strongly dominate over other matrix elements then the 
proportionality is not affected by the many-body effects. 
Indeed, if other weak matrix elements are neglected then 
any many-body expression is a sum of terms with one 
s i/2 — P1/2 weak matrix element in each term. Since all 
of them are proportional with the same proportionality 
coefficient R the proportionality would hold for the sum 
as well. 

Table U shows the ratios of the (6s!/2|-ffw|6pi/2) ma- 
trix elements for cesium atom for different P-odd and 
CP-odd operators of weak interaction ffyv in the rel- 
ativistic Hartree-Fock (RHF) approximation and with 
dominating many-body effects included. The ratios are 
compared with the analytical result presented in the ap- 
pendix. The many-body effects include core polarization 
(CP) and Brueckner-type correlations. Core polarization 
can be understood as the change of the self-consistent 
core potential due to the effect of an external field (weak 
interaction in our case) . Its inclusion above the Hartree- 
Fock level is often called the random-phase approxima- 
tion (RPA). Brueckner-type correlations are the correla- 
tions which can be reduced to a redefinition of the single- 
electron orbitals, replacing the Hartree-Fock ones by the 
Brueckner orbitals (BO) [24]. In both cases of the P-odd 
and CP-odd weak interactions the ratios of the matrix 
elements are stable within the 1-2% accuracy while the 
matrix elements change by up to two times. 



III. PARITY NON-CONSERVATION 

Hamiltonian describing parity-nonconserving electron- 
nuclear interaction can be written as a sum of spin- 
independent (SI) and spin-dependent (SD) parts (we use 
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atomic units: h - 



1): 



_ Gf / Qw 
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(4) 



where Gf ~ 2.2225 x 10 14 a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 

I J ) and 75 = ( _J ~l) are Dirac matri- 

ces, I is the nuclear spin, and p(r) is the nuclear den- 
sity normalized to 1. The strength of the spin-dependent 
PNC interaction is proportional to the dimensionless con- 
stant h which is to be found from the measurements. 
There are three major contributions to h arising from 
(i) electromagnetic interaction of atomic electrons with 
nuclear anapole moment [251 , (i i) electron-nucleus spin- 
dependent weak interaction [26[ , and (iii) combined ef- 
fect of the spin-independent weak interaction and the 
magnetic hyperfine interaction [27| (see, also review [l|). 
In this work we do not distinguish between different con- 
tributions to k and present the results in terms of total 
k which is the sum of all possible contributions. 

Within the standard model the weak nuclear charge 
Qw is given by [11] 
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-0.98777V + 0.0716Z. 



(5) 



Here N is the number of neutrons, Z is the number of 
protons. 

The PNC amplitude of an electric dipole transition 
between states of the same parity \i) = \JiFiMi) and 
|/) = \JfF f Mf) is equal to: 
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where d = ~eJ2i r i 1S the electric dipole operator, and 
F = I + J is the total angular momentum. 

Applying the Wigner-Eckart theorem we can express 
the amplitudes via reduced matrix elements 
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Detailed expressions for the reduced matrix elements of 
the SI and SD PNC amplitudes can be found e.g. in 
Refs. [29] and [30]. For the SI amplitude we have 
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It is convenient to present the amplitude in a compact 
form 
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where C = C(Ff, Jf, Fi, Jj) is the angular coefficient and 
sums Si and S2 do not depend on Ff and F^. 
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where b is the electron part of the SD weak interaction 

(11) 



b = ^Lap(r)>c. 
v2 



Like in the spin- independent PNC amplitude ©, it 
is convenient to present the SD amplitude in a compact 
form 



(J f ,F f \\d S v\\J i ,F i )=J2*S' i . 



(12) 



Here Ci = c(Ff,J n ,Fi) (i = 1,2,3,4) are angular coeffi- 
cients which can be extracted from (fTU|) . Sums S- do not 
depend on Ff and F^. 
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Equation (|T2"j) has more terms than © due to the elec- 
tron vector nature of the nuclear-spin-dependent opera- 
tor. 
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The total PNC amplitude can be presented in a form 
convenient for extraction of the values of x from the PNC 
measurements (see, e.g. |3ll l32l|) 



PNC in cesium 
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= E(1 + R). 



(13) 



Here E is the SI PNC amplitude given by flSJ) and R is the 
ratio of the SD to SI PNC amplitudes. In this work we 
are mostly interested in the values of R. Using compact 
expressions ((9]) and (fl~2|) one can write 
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According to the discussion of the previous section, ratios 
of the SD and SI weak matrix elements do not depend on 
the principal quantum number n. Therefore, using the 
ratios 
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we can get rid of the sums S[, S 2 , S3, S 4 , involving the 
SD matrix elements and replace them with the SI sums 
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The only parameter in (|14[) which comes from numerical 
calculations is the ratio Si/S 2 of two different contri- 
butions to the SI PNC amplitude (see Eqs.flEEJ)). All 
other parameters are given by analytical expressions. 
The ci, c 2 , C3, c 4 and C parameters are just angular coef- 
ficients. The ratios r 1; r 2 , r 3 , r 4 can also be approximated 
by analytical expressions which will be discussed below. 

The expression (|14p can be further simplified in an 
important case of an ns — n's transition (e.g. the 6s-7s 
transition in Cs) . On a few per cent level of accuracy the 
s — p 3 / 2 matrix elements of the SD PNC interaction can 
be neglected (32[. This means that r 2 = r 4 = 0. We also 
have ri — r 3 — r and (|16[) is reduced to 
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Substituting r from (|A10I) leads to 

i? = 4.90(1 - 0.073ZV) Cl5l/ ^ 2 
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Qw is the weak nuclear charge. As in ([TO the only 
parameter which comes from numerical calculations is 
the ratio Si/S 2 . On the other hand, knowing the value of 
R for at least two hfs components of the PNC transition 
is sufficient for extraction of x from the measurements. 



Experimental values for the two different hfs compo- 
nents of the parity non-conserving 6s — 7s transition in 
cesium are 0] 

£pnc(6sf= 4 - 7s_f =3 ) = 1.6349(80) mV/cm, 
£pnc(6s f=3 - 7sir =4 ) = 1.5576(77) mV/cm. (19) 

To extract x we use (fT3"]) and calculate R for these 
two transitions using ([18]) . To do so we note that 
J f = ^ = J n = 1/2, Q w =-73. 19 (see Eq.©) and take 
Si/S 2 = -0.3459 from Ref. (33[. This leads to the system 
of equations 



E(l + 0.06739x) 
E(l - 0.05937^) 



1.6349(80), 
1.5576(77). 



(20) 



The solution for h is >c—0. 382(56). This result is in good 
agreement with the values x=Q. 393(56) from Ref. [Ill ] 
and >f=0. 375(56) from Ref. [12] 1 . Accurate calculations 
similar to what is reported in our previous work for Ba + , 
Yb + and Ra+ [32j lead to the value x = 0.376 which 
is in perfect agreement with the all-order calculations of 
Ref. 121. This is an illustration of the accuracy of the 
analysis based on the ratio of the matrix elements. The 
value x-=0.382 coming from this analysis differs by less 
than 2% from the value x=0.376 coming from the accu- 
rate calculations. This difference is due to two simplifi- 
cations: (a) neglecting the s — p 3 / 2 matrix elements of 
the SD weak interaction, and (b) assuming that the ratio 
of the matrix elements is the same for all single-electron 
states. 

Recent relativistic coupled-cluster calculations of the 
nuclear spin-dependent PNC in Cs |Io[ report the val- 
ues of the SD PNC matrix elements which are about 
30% smaller than those of the all-order calculations of 
Ref. [l2l ]. This is in disagreement not only with this 
work but with all earlier calculations of the SD PNC in 
cesium Given the proportionality of the matrix 

elements of the SI and SD weak interactions discussed 
above, the results of Ref. [10( are also in disagreement 
with all most accurate calculations of the SI PNC in ce- 
sium (see, e.g. 0,0, HI])- The latter calculations have 
accuracy better than 0.5% and are used to test the Stan- 
dard Model in Cs PNC experiment [2| where the accu- 
racy is 0.35 %. 



B. Anapole moment of thallium 

The value of the nuclear anapole moment of thallium, 
extracted from the measurements of the PNC in the 



1 Refs.ri"H and [T3l use different definition of The conversion 
factors are (/ + 1/2) /(/ + 1) for Ref. [JJ] and / for Ref. [JJ]. / 
is nuclear spin, / = 7/2 for 133 Cs. 
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TABLE II: Parameters £(F, F') of the nuclear spin-dependent 
PNC amplitude in the 6pi/2(F) — (>P3/2{F') transitions in 
thallium. 



F 


F' 


this work 


Ref. [16] 





1 


0.947 


1.10 


1 


1 


-0.325 


-0.462 


1 


2 


-0.325 


-0.348 



6P1/2 ~ 6p 3 / 2 transition [151 ] is in disagreement with 
the results of nuclear calculations (see, e.g. ref. fl7j). 
The analysis of the experimental data based on sim- 
ple single-electron approximations gives the value x a = 
—0.22 ± 0.30 [15| . The analysis based on sophisticated 
many-body calculations gives very close value of >t a = 
-0.26 ± 0.27 [lj]. On the other hand, the "best value" 
obtained in nuclear calculations is xr Q — 0.24 17J. To 
extract the sign of >c a from the experiment one needs the 
relative sign of SI and SD amplitudes. Here we show 
that simple analysis with the use of the analytical ratio 
of the matrix elements of the weak interactions supports 
the findings of Refs.[l5|, [l6[ leaving the problem of sign 
disagreement unsolved. 

If we keep only s\/2 — Pi/ 2 matrix elements of the weak 
interaction for both SI and SD interactions, then the gen- 
eral expression (ITBl can be reduced to 



R(F, F') = r 



c(F,F') 
C(F,F')' 



(21) 



where r is one of the ratios (|15l) which corresponds 
to the S1/2 — P1/2 weak matrix elements, c(F,F') and 
C(F, F 1 ) are corresponding angular coefficients. Com- 
paring (fl~3f to the parameterization used in [16] (1Z = 
Im(£l PNC /Ml)) 



K(F, F') = C(Z) [Q w - 6x£(F, F')) , 



(22) 



and substituting (|A10|) into (j21~j) we get for the parame- 
ters £ 



Z(F,F') 



-0.817 (1 - 0.073Z 2 a 2 ) C ^ F F p, y (23) 



Corresponding values of £ are compared in Table ITTI with 
the results of the many-body calculations of [16j ■ We see 
that both calculations give very close results leaving no 
room for a sign error. 



IV. EDM OF ATOMS AND MOLECULES 



Table ITTTl shows the ratio of the S1/2 — P1/2 matrix ele- 
ments of the electron EDM operator (IA15[) to the scalar- 
pseudoscalar CP-odd operator (|A13|) for eight different 
atoms calculated using analytical formulas presented in 
the appendix. The results for cesium and thallium are 
in good agreement with the many-body calculations of 



TABLE III: The ratio of the S1/2 — P1/2 matrix elements 
of the electron EDM operator (|A15|) to that of the scalar- 
pseudoscalar CP-odd operator (|A13[) . Numerical results for 
EDM of Cs, Tl and YbF are also given for comparison. Units: 
d e /(C SP 10~ 18 e cm). (For other isotopes R' = A'R/A) 
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Atom 


Ratios 

Analytical 


Numerical 


37 


8b Rb 


228 




55 


133 Cs 


158 


163° 


56 


138 Ba+ 


152 




70 


173 Yb + 


114 


115 6 


81 


205 rpj 


89 


83 c 


82 


208 Pb 


88 




87 


211 Fr 


83 




90 


232 


75 





"Cs atom, Ref. @. 
6 YbF molecule, Ref. 
C T1 atom, Ref. @. 



Ref. 0|, the result for Yb + is in excellent agreement 
with the many-body calculations of Ref. [36[ for the YbF 
molecule. We stress that the Table compares the ratios 
of the single-electron matrix elements obtained with a 
simple analytical formula to the ratios of the EDMs ob- 
tained with sophisticated many-body calculations. These 
results provide an unambiguous link between the sign 
and value of two different contributions to the EDM of 
atoms and molecules which have a heavy atom from the 
Table HH 



A. EDM of polar molecules 

Polar molecules have strong inter-atomic electric field 
which enhance the effect of electron EDM and lead to 
molecular EDM which are several orders of magnitude 
larger than those in atomic systems. The experimen- 
tal search is in progress for YbF [18], PbO [2H l22|. and 
ThO [23| while other molecules are also discussed in the 
literature (see, e.g. [34]). Interpretation of the measure- 
ments requires molecular calculations. Table IIVI shows 
the results of most recent or most accurate calculations 
for some polar molecules. More detailed data are pre- 
sented for YbF molecule for which the EDM measure- 
ments were recently reported [l8j]. The effects of elec- 
tron EDM and scalar-pseudoscalar CP-odd interaction 
are considered. The results are presented in terms of the 
CP-odd parameters Wd and W c 



W d = (* \H e \* )/d e , 



Wr. = (tfnlff 



SP 



SP 



(24) 
(25) 



To compare with other works one should keep in mind 
that most of them present Ws instead of W c , where 
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(26) 
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TABLE IV: CP-odd interaction constants W d , and W c for 
some polar molecules and their ratios. Effective electric 
field E a ff is presented together with W d ({H e ) = W d d e — 
— E e fjd e ). The results of present paper are shown in bold. 







w d 


W c 


w d /w c 


Molecule 


(GV/cm) 


(10 24 Hz/e cm) 


(kHz) 


(10 18 /e cm; 


BaF 


6.1" 


-1.5" 


-10 


152 




6.1 


-1.5 


-10 6 


152 


YbF 


31 c 


-7.5 C 


-59 c 


127 




19 d 


-4.6 d 


-41 d 


112 




26 e 


-6.3 C 


-55 


114 




26 f 


-6.2 / 


-54: f 


115 




24 s 


-5.8 9 


-54 9 


108 




25 h 


-6.1 ft 


-53 


114 




32" 


-7.7 a 


-68 


114 




21 


-5.2 


-46 6 


114 


HgF 


95" 


-23° 


-226 


90 


PbF 


-31 a 


7.5° 


85 


88 


PbO o(l) 


25' 


-6.P 


-69 


88 


PbO 5(1) 


33 1 


-8.0 s 


-91 


88 


PbO o(l) 


23° 


-5.6 a 


-64 


88 


ThO 


104" 


-25° 


-336 


75 


ThF+ 


90" 


-22 a 


-290 


75 



"Meyer and Bohn, Ref. Eg. 
6 Nayak and Chaudhuri, Ref. (H). 
c Kozlov and Ezhov, Ref. [jjj. 
d Titov et al, Ref. 
e Kozlov, Ref. [39| . 
^Quiney et al, Ref. |36H . 
fParpia, Ref. [2(J. 
h Mosyagin ai, Ref. [H. 
"Petrov et al, Ref. [3. 

The constants k s and C SP of the strength of the CP-odd 
interaction are related by Zk s = AC SP . Factor 2 in the 
definition of Ws ([26]) is absent in some of the papers. 

Most of the calculations of the molecular EDM in- 
clude only one of the CP-odd effects: that of the elec- 
tron EDM or the scalar-pseudoscalar interaction. We 
use the relations between the matrix elements of the two 
CP-odd Hamiltonians to fill the gaps in the table. Corre- 
sponding results are shown in bold. For example, accord- 
ing to Ref. [H the effect of electron EDM in the YbF 
molecule is Wd = — 15 x 10 24 Hz/e cm. Using the ratio 
W d /W c = 114 x 10 18 /e cm for Yb+ from Table [[n] we 
found that the effect of the scalar-pseudoscalar interac- 
tion is W c = -132 kHz. 



B. Extraction of the electron EDM and the 
parameter of the scalar-pseudoscalar CP-odd 
interaction from the EDM measurements for YbF 
and Tl 

Recent measurement of T and P violation in YbF 
molecule [l8| combined with the Tl EDM measure- 
ment [l9| and the data from Table IIIII allows to obtain 
independent limits on the electron EDM and the CP- 



violating interaction. The values of the electron EDM 
extracted from the experimental data for Tl 19] and 
YbF [l8| under assumption that there is no other contri- 
bution to atomic/molecular EDM read 

Tl : d e = (6.9 ± 7.4) x 1(T 28 e cm, (27) 
YbF: d e = (-2.4 ±5.7 ±1.5) x 10~ 28 e cm. (28) 

In fact, there are other contributions from the CP-odd 
electron- nucleus interactions. Here we consider only the 
scalar-pseudoscalar interaction (lAllj) . Other contribu- 
tions should be small due to the constrains obtained from 
the EDM measurements for mercury (43| . Using data 
from Table IIIII one can rewrite (|27I28[) as 

Tl : d e + 1.2 x 1(T 20 C SP = 

(6.9 ± 7.4) x 1(T 28 e cm, (29) 
YbF : d e + 8.8 x 1CT 21 C* SP = 

(-2.4 ± 5.7 ± 1.5) x 1(T 28 e cm. (30) 

Solving these equations for d e and C leads to 

d e = (-2.8 ± 3.0) x 10" 27 e cm, (31) 
C SP = (3.0 ±3.0) x 10~ 7 . (32) 

The limit pip for d e is significantly weaker than those 
presented in (|2T|) and l[2"8]). This is because the values 
of the ratios of the matrix elements of the two operators 
are very close for Tl and Yb (see Table HIT)) . The most 
fortunate case leading to strongest limits for d e and C SP 
would correspond to very different values of the ratio, 
preferably with the different sign. However, the formulas 
(IA15|) and (|A13j) and the data in Table |HT] show that this 
is not possible. The ratio is always positive and slowly 
decrease with Z. A slightly different value of the ratio 
should be expected for the EDM of radium in the excited 
metastable 7s6<i 3 D2 state. The EDM of radium in this 
state is strongly enhanced [44| due to proximity of the 
7s7p 3 P° state (the energy interval is ~ 5 cm -1 ) and 
is proportional to the (7p3/2\Hcp\6d 3 / 2 ) matrix element 
of the CP-odd interaction Hep- However, the value of 
this matrix element is strongly dominated by the core 
polarization effect which in turn is mostly due to the 
s i/2 — P1/2 matrix elements between core and excited 
states. Therefore, the ratio of the two contributions is 
not very different from what was considered above for Tl 
and YbF. 



C. EDM of Tl and Fr 

EDM of thallium due to electron EDM and scalar- 
pseudoscalar CP-odd interaction was recently calculated 
in our paper [^]. The ratio of the two contributions 
was found to be 83d e /(C SP 10~ 18 e cm) which is in 
good agreement with the analytical result of this pa- 
per: W d /W c = 89/(10" 18 e cm). The EDM enhance- 
ment factor is —582 0] which is in very good agree- 
ment with the value —585 found in earlier calculations 
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by Liu and Kelly @. The EDM enhancement factor 
for thallium was also calculated in more recent work of 
Ref. and found to be -466 (g?ti = — 466d e ) which is 
about 25% smaller than the results mentioned above. 
The EDM of thallium due to the scalar-pseudoscalar 
CP-odd interaction calculated bythe same group earlier 
(rfxi = -4.06 x 10~ 18 C SP e cm pj) is also smaller than 
in other calculations of Ref. Q and (20j . On the other 
hand, the ratio of the two contributions is 115 (in units 
of d e / (C SP 10~ 18 e cm)) which is significantly larger than 
the analytical and numerical values of 89 and 83 (same 
units, see Tabic lUT)) . 

EDM enhancement factor for francium was calculated 
in Ref. [i| and found to be 910(46): rf(Fr) = 910(46)4. 
This value is in good agreement with the value of 894.93 
found in more recent calculations [H(J. Using the value 
d(Fr) = 910(46)4 and the ratio of the electron-EDM ma- 
trix element to the scalar-pseudoscalar interaction pre- 
sented in Table Mil we can now reconstruct another con- 
tribution to the EDM of francium: 

<2(Fr) = 11 x 10- 18 C SP e cm. (33) 

It is interesting to note that the calculation of the nu- 
clear spin-dependent PNC amplitude between hypcrfinc 
components of the ground state of francium presented 
in Ref. p9l | involves the same sum as the expression 
@ for the EDM of the atom. Therefore, we can re- 
construct the EDM enhancement factor for Fr using the 
results of Ref. [2{| and proportionality between matrix 
elements of the nuclear spin-dependent interaction (|A9j) 
and the electron EDM operator (|A15I) . This leads to 
d(Fr) = 854d e which is in good agreement with the value 
d(Fr) = 910(46)d e of Ref. [4§|. This is a good consistency 
test of both calculations. 

Reading the papers citing our works (see, e.g. [(j l50j) 
reveals the need to clear some points about our method 
of calculations. For example, the atomic electric field in- 
teracting with electron EDMs is calculated in [H and [49[ 
as a derivative of the total potential which includes both 
nuclear and electron parts. The formula E = Zer/r 3 
for the leading contribution to the atomic electric field 
presented on first page of [49j may make an impression 
that only nuclear field is included. However, the formula 
(2) few lines below clearly includes screening functions 
Q(r) and P(r) for both the nuclear Coulomb and the 
external electric field. By the way, the inclusion of the 
electron electric field change the matrix elements of the 
electron EDM for thallium by 0.4% only. This is because 
main contribution comes from short distances where the 
electron electric field is small since the electron potential 
rapidly tends to a constant inside the Is orbital. There- 
fore, we include only nuclear electric field in the ana- 
lytical analysis (formula (|A15|) in the appendix) while 
keeping both contributions in the numerical calculations. 

The authors of Q claim that atomic core is strongly 
contracted in the V N ~ 3 starting approximation used in 
our calculations @. In fact, it is not. FigQ]shows the out- 
ermost 5g?5/2 core function of Tl calculated in the V N ~ 3 




1 2 3 4 5 

R (a.u.) 



FIG. 1: The 5d 5 /2 core function of Tl calculated in V 
(solid line) and V N_1 (dashed line) approximations. 



and V approximations. The difference between the 
functions is very small. This is due to the fact that the 
valence 6s and 6p electrons are located outside of atomic 
core. Their charge distribution creates almost constant 
potential and no electric field inside the sphere where all 
inner electrons are located. Therefore, the valence elec- 
trons have practically no effect on the core wave func- 
tions (see [45[ for a detailed discussion). The change is 
even smaller for other core functions. The core functions 
enter the configuration interaction (CI) Hamiltonian via 
core potential V core to which all core electrons contribute 
(from Is to 5d). The difference for V core in the V N ~ 3 
and V N ~~ l approximations is very small [45| . Moreover, 
the corresponding corrections to the configuration inter- 
action (CI) Hamiltonian have been included in Q using 
the many-body perturbation theory approach. 

As it is well known, the eigenstates of a Hamiltonian do 
not depend on the basis one uses. The valence states are 
indeed different in the V N ~ 3 and V N ~ 1 approximations. 
However, this should have no effect on the final results 
as long as the saturation of the basis for valence states is 
achieved. There are only two conditions the basis states 
must satisfy: (a) they must be orthogonal to the core, 
and (b) they must constitute a complete set of states. 
Therefore, it does not matter whether valence states are 
calculated in the V N ~ 3 or V^ -1 potential or by any other 
method (e.g., a popular basis ip n (r) = r n ip (r) [46j), the 
final results should be the same. If there is any differ- 
ence in the results, the most likely reason for this is the 
incompleteness of the basis set. 

In spite of no difference in final results there is a good 
reason for the use of the V N ~ 3 approximation - the sim- 
plicity and good convergence of the many-body pertur- 
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bation theory (MBPT) for the core-valence correlations. 
This is the only approximation which has no subtraction 
diagrams [13]. To be more precise, we should state that 
the condition for the absence of the subtraction diagrams 
is that the potential used to calculate single-electron core 
states and the potential of the core in the CI Hamiltonian 
are exactly the same. In the case of thallium atom treated 
as a three-valence-electrons system this corresponds to 
the V N ~ 3 approximation. When any other approxima- 
tion is used one has to include the subtraction diagrams. 

Large energy denominators suppress the value of the 
correlation terms in the V N ~ 3 approximation ensuring 
good convergence of the MBPT [45| . There must be large 
cancelation between subtraction and other diagrams to 
ensure the same final results if any other initial approx- 
imation is used. This is very similar to the well known 
fact that the Hartree-Fock basis is the best choice for any 
MBPT calculations. Zero-order results might be better 
in some other approximation, however, strong cancela- 
tion between subtraction and other diagrams would lead 
to poor convergence of the MBPT. 

The authors of Q claim that the major drawback of 
our work is the inclusion of the core- valence correla- 
tions in the second order only. However, the correlations 
between the valence electrons and core electrons below 
the 6s state are small which is evident from the fact that 
their inclusion change the EDM of Tl by 3% only @. 
Therefore, only the correlations between three valence 
electrons should be treated to all orders. This is done 
in Q to a very high precision using the CI technique. 
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Appendix A: Matrix elements of the weak 
interaction 

1. Wave function 

We use single-electron wave functions in a form 
1 / f v (r)Q(n)ji m 



ip(r) 



njlm 



r \ iag v (r)Cl(n)ji„ 



(Al) 



where n is the principal quantum number and an index 
v replaces the three-number set n, j, I; a is the fine struc- 
ture constant. 

We need S\/2 and p\/2 wave functions inside the nu- 
cleus to calculate matrix elements. Following Ref. [5l[ 
we assume uniform distribution of electric charge inside 
a sphere of radius Rn- Taking formulas for the S1/2 and 
P1/2 wave functions from Ref. [5l| and keeping terms up 
to Z 2 a 2 we come to the following expressions 



(A2) 



CONCLUSION 



f s (x)=A s x \l--Z z a'x z y 
/pi/aO) = A p Za 2 — 1 - — I , 



9pi/A x ) = A p x ( 1 - cZ 2 a 2 x 2 



We have demonstrated the proportionality rela- 
tions between the nuclear-spin dependent and spin- 
independent PNC effects in atoms and the scalar- 
pseudoscalar and the electron EDM contributions to the 
EDM of atoms and molecules. The relations are given 
by the simple analytical formulas and can be used to 
express one symmetry breaking effect through another. 
Using these relations and accurate calculations of the 
spin-independent PNC we have confirmed earlier inter- 
pretations of the nuclear anapole measurements in cesium 
and thallium. We have also confirmed the ratio of the 
scalar-pseudoscalar CP-odd and electron EDM contribu- 
tions to the EDM of Cs, Tl and other atoms and some 
polar molecules. Using the relations we found the scalar- 
pseudoscalar contribution to the EDM of francium atom 
and the scalar-pseudoscalar and electron EDM contribu- 
tions to the EDM of many polar molecules. Using exper- 
imental limits on EDMs of thallium and YbF we found 
model-independent limits on the electron EDM and the 
constant of the scalar-pseudoscalar CP-odd interaction. 



Here x = r/i?jv, Rn is nuclear radius. A s and A p are 
normalization factors. They can be found by matching 
(|A2|) to atomic wave functions at short distances outside 
of the nucleus (see, e.g. [Hj]). The use of the semiclassical 
wave functions presented in [52| leads to approximate 
expressions 



A s = 
A n = 



(1 - 0.375Z 2 a 2 )V2T(2 7 + l)v\* 

VZ(2ZR N )~> 
(1 - 0.375Z 2 a 2 )r(2 7 + l)j/i- 5 ' 



(A3) 



where 7 = y/l — Z 2 a 2 , T is gamma-function, v is effec- 
tive principal quantum number: single-electron energy of 
state n is e n = —l/(2is 2 ). Note that normalization fac- 
tors (IA3|) are needed only to compare the matrix elements 
of the electron EDM to other matrix elements. For cases 
not involving electron EDM there is exact cancellation of 
the normalization factors. Therefore, their uncertainty 
does not contribute to the uncertainty of final results. 



9 



2. PNC matrix elements 



3. Scalar-pseudoscalar CP-odd interaction 



Matrix element of the SI PNC Hgi interaction (first 
term in Q) is 



/ it i \ -GfQw 

/ (/i52 - 9ih)p{r)dr. 



(A4) 



Substituting (|A2j) for the Sx/2 ~ Pi/2 transition and as- 
suming uniform nuclear charge distribution we get 



/ itt i \ -GfQw 
(si/al-Hsibi/a) = x 

ap A s A p R N {1 _ QMZ 
3 



(A5) 



2 a 2 ) 



Here A s and A p are given by (|A3[> and 

3 

Po = 



Matrix element of the SD PNC b operator (fTTj) is 

(ki||6||k 2 ) = — i — 7^a(— /ci||Ci||k2) x 
V2 



(A7) 



((«2 - «i + l)/i52 - («x - K 2 + 1)31/2) p(r)dr. 



The reduced matrix element of the spherical harmonic 
C fc is 



(«al|C k ||« 6 ) = (-iy» +1 /V(2ja + l)(2j6 + l) 

Xe(ia + J6 + *)(_{/2 1/2 o)' (A8) 

1, if x is even 



0, if x is odd 



For the Sx/ 2 — P1/2 transition and uniform nuclear charge 
distribution we have 



(S1/2II&IIP1/2) = ~i— F -^- x 



V2 



a\/ - Po A s A p R N (1 - 0.41ZW 



(A9) 



The Hamiltonian of the scalar-pseudoscalar electron- 
nucleon (T,P)-odd interaction can be written as [l| 



H sp = i^AC sp WN (r), (All) 

where Gf is the Fermi constant, A = Z + N is the nu- 
clear mass number, Z is the number of protons and N 
is the number of neutrons. C SP = (ZC^ P + NCf p )/A : 
where C° p and Cf p are the parameters of the scalar- 
pseudoscalar (T,P)-odd interaction for protons and neu- 
trons, 7„ are the Dirac matrices. 



(a\H sp \b) = -^faAC s 
v2 



x<5_ 



(A6) Substitution of the functions (|A2I) leads to 



( S l/2|^ SP bl/2) 



(fa9b + g a fb)pNdr. (A12) 



(A13) 



G F AC SP 
3V2 



aA s A p poR N (1 - 0.56ZV 



Here A s and A p are given by (|A3[) and po is given by 



4. Electron EDM 

The Hamiltonian for the electron EDM interacting 
with internal atomic electric field E int can be written 
as 52] 



H e = -4(70 - 1)S - Eint, 



(A14) 



where 



S = 



a 
cr 



The ratio of (1A9I) to (IA5I) is 



and Ej nt is internal atomic electric field. The Sx/2 — P1/2 
matrix clement can be found in Rcf . 53] : 



(Sl/2lNPl/2>/(Sl/ 2 |ffsi|Pl/2> « (A10) 



( ^5F^ (1 - - 073zV) 



(si/ 2 |i?ebl/2 



4Z 3 a 2 d, 



7(472 - l){v s v p y/* 



(A15) 
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